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We study the magnetism and lattice distortion of Fei+yTei_xSex with a general field theory formalism, 
motivated by recent neutron scattering experiments. Besides the Ising nematic order parameter which can be 
naturally defined in the 1111 and 122 materials. We show that the coUinear order observed in Fei+yTe materials 
is stabilized by an extra Ising order parameter which doubles the generalized unit cell of the system. By tuning 
this Ising order parameter, one can drive a transition between collinear and noncoUinear spin density wave 
orders, both of which are observed in the FeTe family. The nature of the quantum and classical phase transitions 
are also studied in this field theory formalism, with lattice strain tensor fluctuation considered. 



The Iron-superconductor, for its potential to shed new light 
on the non-BCS type of superconductors, has attracted enor- 
mous interests since early last year So far most studies have 
been focused on the 1111 materials and 122 materials. The 
1111 materials have chemical formula LnFcAsO, and Ln 
represents one of the lanthanides; The 122 materials have 
chemical formula MFe2As2, and M usually represents one of 
the alkaline earths. At low temperature, it is confirmed that 
both 1111 and 122 materials develop (tt, 0) spin density wave 
(SDW) nil Hi, Ever since its discovery, the SDW 

in 1111 and 122 materials raised a debate about whether this 
phenomenon should be described by an itinerant electron pic- 
ture with nesting fermi surface, or a local moment Heisenberg 
model. This debate is partly due to the fact that the order- 
ing wave vector (tt, 0) connects two pieces of nearly parallel 
fermi surfaces, and relatively small magnetic order moment 
of 111 1 materials. However, in the relatively less extensively 
studied 1 1 material FcTc family, the answer seems to be more 
clear, as the order wave vector of these materials is exactly or 
close to (7r/2, 7r/2) HI, V2\, which is far away from the nest- 
ing wave vector according to LDA calculation 122 1. Also, the 
magnetic order moment can be as large as ll ill , which 
can be naturally understood in terms of local moment picture. 

In the low temperature phase of Fci+yTc, besides the spin 
density wave (SDW) with wave vector (7r/2,7r/2), a mono- 
clinic lattice distortion ll ill has also been identified. By in- 
creasing the excess of Fe in this system, the magnetic order 
becomes incommensurate lfl2ll . with order wave vector devi- 
ates from but still close to (7r/2,7r/2). Also, by replacing 
Te with Se, both SDW and lattice distortion are suppressed. 
In Ref. 10], the authors used a local moment Heisenberg 
model with the nearest, 2nd nearest and 3rd nearest neighbor 
( Ji — J2 — J3) interactions on the square lattice to describe 
this system. This simple model captures the correct phases 
observed experimentally and is also supported by results from 
LDA calculations llTol] . But a model independent general field 
theory formalism is also necessary, in order to for instance ex- 
plore the physics beyond the Heisenberg model, and to study 
the nature of the phase transitions. In this work we will sys- 
tematically study the Ginzburg-Landau field theory for the 



magnetism and lattice structures of FeTe compound and its 
relative Fei+yTei_xSex. 

A general Ginzburg-Landau-Hertz-Millis theory for the 
1111 and 122 materials was studied in Ref. IiITIi. M otivated 
by experimental facts, it was proposed in Ref. 11141 [Till that 
the true magnetic order in 1111 materials (say LaFeAsO) is 
developed through two separate steps: first an anisotropic an- 
tiferromagnetic (AF) correlation is developed at a relatively 
higher temperature, and then the long range spin density 
wave (SDW) order is developed at lower temperature. The 
anisotropic AF correlation described by an Ising order param- 
eter breaks the reflection symmetry of the lattice along the 
axis X = y, and leads to a tetragonal-orthorhombic lattice dis- 
tortion. More specifically, this Ising order parameter can be 
defined as $ ~ 01 • (/)2, 0i and (j>2 are Neel order parameters 
on two sublattices of the square lattice. 

In the FcTc family, the SDW and lattice distortion are more 
complicated, and presumably the anisotropic AF correlations 
can also develop before the true long range SDW. The SDW 
order pattern contains four Neel orders on four sublattices 
I nil (Fig- [D; 0a with a = 1 — 4. Two Ising order param- 
eters which obviously break the lattice reflection symmetry 
are $1 ~ 01 • 03 and (f>2 ~ 02 • 04- 'I'l and $2 can develop 
long range order before 0a does, and the symmetry of en- 
ables them to distort the lattice from tetragonal to monoclinic 
once they develop long range order. However, $1 and $2 are 
not independent Ising order parameters i.e. in the free energy 
there should be a term which breaks the degeneracy between 
<f>i<l'2 > and <f>i$2 < 0. Later on we will show in the 
ground state $i<f>2 < 0, and the energy barrier between states 
with (f>i$2 < and $i(f>2 > is estimated to be about Jf / J3 
for the Ji ~ J2 — J3 model. 

Let us first study the structure transitions, and tentatively 
ignore the magnetic orders. In the high temperature symmet- 
ric phase, the Fe-Te plane enjoys the following symmetries: 
Pz (8) Tx, Pz iXi Ty, Pz (8) Fx, Pz Py Pz and Fx are reflec- 
tions z ^ —z and i; — > — i with origin located on one of the 
Fe atoms, Tx and Ty are translations. Although the true unit 
cell of the Fe-Te plane always involves two Fe atoms, one can 
define a one Fe unit cell based on the generalized translations 
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FIG. 1: The collinear magnetic order of FeTe family. The four sublat- 
tices form Neel orders separately, and the thick lines represent links 
with stronger AF correlations, which doubles the unit cell defined 
by generalized translation operators T'^ and Ty. It is possible that 
the anisotropic AF correlations develop long range order before the 
actual SDW pattern. 



Tx' = Pz Tx and Ty' = (g) Ty, and we will use this 
generalized unit cell hereafter. The lattice distortion observed 
in experiments lull breaks all of the symmetries including T^' 



and Ty . The lattice distortion favors the antiferromagnetic 
correlations along the zig-zag stripes, and there are in total 4 
different degenerate zig-zag configurations. We can describe 
these zig-zag AF correlation stripes with two Ising order pa- 
rameters (Ji and (72. CTi = ±1 represents two degenerate zig- 
zag correlation stripe patterns along the x = y direction, and 
(72 = ±1 represents two zig-zag correlation stripes along the 
X = —y direction. One can write down a Ginzburg-Landau 
free energy for the coarse-grained modes of (Ta according to 
the symmetry of the system: 
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Here we need U2 > 2uo to make sure that in the ordered phase 
(o'i)(o'2) = 0, which is the case in the real system. 

Another more intuitive way to describe the (Ta variables, is 
to introduce complex field 'i' = ai + i(J2, and the free energy 
Eq. [T]can be rewritten as 



r\^\' + + gii'f + * ) + 0(*^).(2) 



The free energy up to the forth order of aa and ^ enjoys an en- 
larged Z4 symmetry, which can be viewed as an Z4 anisotropy 
of an XY or U(l) quantity 5'. It is well-known that the Z4 
anisotropy at the 3D XY transition is irrelevant 1IT9I1 . so our 
first conclusion is that, in the most simplified situation the fi- 
nite temperature transition of era is a 3D XY transition. In Eq. 
[Hand Eq. |2]The terms at the 6th order or higher will break 
this Z4 symmetry, however, those terms will also be irrelevant 
at the 3D XY transition, so we will ignore them from now on. 



According to the symmetry, the Ising fields cti and (T2 cou- 
ple to the Ising field ^^^f&i^— $2in the following manner: 



(3) 



The Ising field $ describes the anisotropic AF correlation 
which breaks no translational symmetry, but breaks the reflec- 
tion symmetries Pz ® Px and Pz ® Py. The full free energy 
reads: 



r2$' 



(4) 



The ellipses include all the quartic terms. Both ri and r2 are 
tuned by temperature. By minimizing Eq. |4] one can obtain 
the mean field global phase diagram plotted against r = ri + 
r2 and Ar = — r2 (Fig. |2]2). There are in total three 
different regions: 1), when ri 3> ^2, there are two transitions, 
with one Ising transition for $ and another Ising transition 
for one of the aa, since in this region $ has a much stronger 
tendency to order compared with 0-^; 2), when ri ~ r2, there 
is an intermediate region with first order transition, and $ and 
(Ta will order at the same temperature. The first order nature of 
this transition is due to the cubic coupling term F$,cr. 3), when 
ri <C 7*2 there is still one single transition, but this transition 
is second order. This transition can be studied by integrating 
out <1> from free energy Eq. |4] and the resultant free energy for 
(Ta is the same as Eq. |2l which describes an 3D XY transition 
as was discussed previously. 

So far we have ignored the lattice elasticity, which will po- 
tentially change the nature of the transitions discussed above 
IitIi . In region 1) discussed in the previous paragraph, the 
Ising field $ will couple to the uniaxial shear strain field of 
the lattice: 



dyUy), 



(5) 



u represents the local displacement vector of the lattice. By 
integrating out the displacement vector u, we can show that 
the fluctuation of the displacement vector u will increase the 
effective spatial dimension at the transition of $, and drive 
the transition of $ mean field like, as was studied carefully in 
Ref. UtIi . At the mean field transition the specific heat curve 
will have a discontinuity instead of a peak, as was observed in 
electron-doped 122 materials Unlike the field in region 
1) the field (Ta that breaks the translational symmetry will only 
couple to the bulk strain field in the following way: 
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We can estimate the scaling dimension of the coupling con- 
stant fa-,u as the following: The scaling dimension of |(Tap 
at the 3D Ising transition is A[|crap] = 3 — l/u, and the 
scaling dimension of displacement vector u is A[u] = 1/2. 
Therefore the dimension of is A[/o.,u] ~ —3/2 + 1/i^ = 
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FIG. 2: The global phase diagrams for Eq. |4]and Eq. 1101 plotted 
against r — ri +r2, and Ar — ri —r2. ForEq. ID Fig. a represents 
the case without considering lattice elasticity. In this case the red line 
is a first order transition; the golden line is a 3D XY Wilson-Fisher 
transition; the green line and blue line are 3D Ising transitions for <I> 
and aa respectively. For Eq. |4] Fig. b represents the phase diagram 
with coupling to the lattice strain tensor, the 3D XY transition and the 
3D Ising transition of aa become first order, the 3D Ising transition 
of $ becomes a mean field transition. For Eq. [TO] Fig. a and Fig. 
b represent the case without and with considering itinerant electron 
particle-hole excitations respectively. If a two dimensional model is 
taken, the nature of the transitions is almost the same as Eq. |4] 



aismg/2 > 0. So thls coupling will lead to a relevant per- 
turbation at the 3D Ising transition of aa- This relevant per- 
turbation will likely drive the transition first order 12111 . This 
argument of scaling dimension is essentially the same as the 
well-known Harris criterion in the disordered system. 

In region 3), as was shown before, after integrating out <!>, 
the transition is described free energy Eq. [21 which describes 
a 3D XY transition. The 3D XY transition is also affected by 
the lattice elasticity. First of all, j^fp is also coupled to the 



bulk strain field ~ l^^l {^xUx + ^yUy + ^^zUz), be- 
cause axy = 2/v ~ D < this coupling will not introduce 
any relevant perturbations. However, ct^ — erf Re[^'^] will 
couple to the shear strain field OxU^ — dyUy. To check whether 
this coupling is relevant or not, we just need to compare the 
scaling dimension of Re ['3/^] to Z)/2 at the three dimensional 
XY transition. At the 3D XY critical point, the scaling dimen- 
sion A[Rc[4'2]] ^ 1.234 < 3/2 HI, hence this coupling will 
indeed generate relevant perturbation at the 3D XY transition, 
which likely drives the transition a first order one. Therefore 
after considering the lattice elasticity, the global phase dia- 



gram is modified to Fig. |2f). 

Now let us consider the SDW order, described by four Neel 
order parameters 0^ with a = 1 — 4, and let us ignore the Ising 
fields (Ja first. Besides the ordinary kinetic terms in the free 
energy, several linear spatial derivative terms are also allowed 
by the symmetry of the system, therefore a most general free 
energy for the Neel order parameter 0^ reads: 
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(7) 



The z direction dispersion has been ignored. This theory can 
also be viewed as a low energy field theory for the lattice Ji — 
J2 — J3 Heisenberg model, with f3 ^ Ji, and g ^ JI/J3. If 
<f>i$2 < 0, An infinitesimal Ji inevitably drives the system 
to an incommensurate order, with ordering wave vector Q = 
{TT/2+6q,T:/2+6q),dq ~ J1/J3. However, if $1 $2 > 0, the 
minima of the dispersion remains atQ^ (7r/2,7r/2) for small 
enough Ji . The energy difference between the two states with 
(f>i$2 < and <f>i$2 > is about J1/J3, or in other words 
there is an effective coupling <f>i$2 in the free energy, which 
is proportional to Jf j J 3. 

In the phase with <i>i$2 < The collinear commensurate 
order is stabilized by coupling the SDW order parameter to 
the Ising fields a a- Again, the symmetry of the system allows 
the following coupling; 



F 1 



7(Ti((^i • 02 + (^2 • '/'S - '/'S • 04 - 01 • 04) 
+ 7O'2(01 ■ 02 - 02 ■ 03 - 03 ■ 04 + 01 ■ 04)- (8) 

After diagonalizing the quadratic part of the entire free energy 
(Eq. [8]l, we can see that if <I>i > and (1>2 < the system 
favors {(J2) ^ while (cri) = 0. Therefore there is indeed an 
effective coupling F^ rr, as Eq. [3] 

Whether the system develops incommensurate or commen- 
surate order, depends on the competition between the terms 
with linear spatial derivatives, and the terms with coupling 
to aa- When ^{aa) dominates Jf/Ja, the commensurate 
collinear order is stabilized, otherwise the system develops the 
incommensurate coplanar order Assuming $ = <i>i > 0, and 
hence at low temperature a2 orders while ai remains disor- 
dered, a transition between incommensurate and commensu- 
rate orders can be driven by tuning a = {a^}- Experimentally 
this transition can be achieved by changing the excess of Fe 
1 1211 . If we diagonalize the quadratic part of the SDW free en- 
ergy, the mode with the lowest minimum takes the following 
form: 

E{q) ^ -gm\ + .h{ql + g^^) - ^ 1^ + {qx - qy? -{9) 

One can check that by tuning a, there is a transition between 
commensurate and incommensurate order, and the incommen- 
surate wave vector is a continuous function of a- The com- 
mensurate collinear SDW order has ground state manifold S*^ 
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times lattice transformation degeneracy, while the incommen- 
surate coplanar SDW order has ground state manifold jZi. 
The transition between the these two phases can be viewed 
as condensation of an U(l) degree of freedom lIlSll . because 
the manifold is locally 5^ x S*^. However, as was shown 
in Ref. ifisll . the condensation of this U(l) degree of free- 
dom does not belong to the XY universality class, because the 
spin-wave fluctuation of the collinear SDW order parameter 
at the transition will again effectively increase the spatial di- 
mension of the XY transition, therefore the collinear/coplanar 
phase transition always belongs to the mean field universality 
class, for both 3D classical and 2H-ld quantum cases jlS I. 

It is also interesting to study the quantum phase transitions 
of a a and $ at zero temperature, and we will take a simple two 
dimensional model to emphasize the physics within the Fe- 
Te plane, and also partly motivated by the fact that the band 
structure calculated by LDA has a much weaker z direction 
dispersion for FeTe and FeSe compared with the 122 com- 
pounds I22II . The Lagrangian describing the quantum phase 
transitions is very similar to the free energy Eq. |4l 
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[TOl The coupling between IctqI and fermion density will lead 
to a term 



+ K +4<?^ +r2)|<I> 

+ /$.a<f(^?-a2) + . 



(10) 



The phase diagram is actually quite similar to Fig. |2l with ?- 
tuned by replacing Te with Se in Fci+yTc, as shown experi- 
mentally illlll2ll . Phase diagram Fig. |2]i corresponds to the 
case without fermi pockets, or phase transitions in a supercon- 
ductor with fully gapped fermi surface, and the nature of the 
phase transitions described by Eq. [TO]is the same as the phase 
transitions in Eq. H 

The fermi surface of the 1 1 materials have been observed 
by ARPES 1E3I1 . and the structure qualitatively agrees with 
the LDA calculation. If the fermi surfaces are not gapped, 
the coupling between the bosonic order parameters a a and $ 
with the particle-hole excitations around the fermi surface will 
change the nature of the quantum phase transitions, or more 
specifically change the dynamical exponent z. Phase diagram 
Fig 12) corresponds to the case with ungapped fermi pockets. 
The green line becomes a z = 3 quantum phase transition, as 
discussed in Ref. |3 13], because the order $ carries zero 
momentum, and the particle-hole excitation around the fermi 
surface will induce a decay term of $: 



(11) 



A similar decay term does not exist at the quadratic level in the 
Lagrangian of o-q, because CTq cjirries momentum (7r/2, 7r/2). 
Based on the LDA calculation I22II and also ARPES measure- 
ments the momentum (7r/2, 7r/2) does not connect two 
pieces of the fermi surfaces. However, the particle-hole exci- 
tation can also modify the quartic terms of the Lagrangian Eq. 



in 



(12) 



which is a relevant perturbation at the 3D Ising transition 
due to the positive aising- On the left side of Fig. |2) 
with Ar < 0, the transition was a 3D XY transition with- 
out fermi pockets. The Fermi pockets will generate a term 
Re[vI/2]^ -MRe[v]/2]_^ which is again relevant at the 3D 
XY transition. Therefore the red line in Fig |2f) is a first order 
transition for a 2 dimensional quantum model. 

In summary, we have shown that Fci+yTc materials are 
characterized by two Ising order parameters: one Ising ne- 
matic order which has been found in 1111 and 122 mate- 
rials and the other Ising order which doubles the unit cell 
of systems. We have also calculated the phase diagram and 
the nature of phase transitions. Unlike long-range SDW or- 
ders which are easily destroyed by replacing Te with Sc in 
FeTei-xScx, the Ising orders may survive in larger doping 
regions llT4l[l6ll and coexist with superconducting states. Ex- 
perimentally, it is possible to detect the Ising orders from their 
broken symmetries, for example, using polarized light in an- 
gle resolved photoemission spectra. 
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